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Abstract 



We study electron transport through a small metallic island in the per- 
turbative regime. Using a diagrammatic real-time technique, we calculate the 
occupation of the island as well as the conductance through the transistor at 
arbitrary temperature and bias voltage in forth order in the tunneling ma- 
trix elements, a process referred to as cotunneling. Our formulation does not 
require the introduction of a cutoff. At resonance we find significant modifica- 
tions of previous theories and quantitative agreement with recent experiments. 
We determine the renormalization of the system parameters and extract the 
arguments of the leading logarithmic terms (which can not be derived from 
usual renormalization group analysis). Furthermore, we perform the low- and 
high-temperature limits. In the former, we find a behavior characteristic for 

the multichannel Kondo model. 
73.40. Gk,72.23.Hk,73.40.Rw 
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I. INTRODUCTION 



Electron transport through small metallic islands is strongly influenced by the charging 
energy associated with low capacitance of the junctions^!. A variety of single-electron 
effects, including Coulomb blockade phenomena and gate-voltage dependent oscillations of 
the conductance, have been observed. If the conductance of the barriers is low 

a = V(47r 2 e 2 i? T ) < 1 (1) 

they can be described within the "orthodox theory" which treats tunneling in lowest order 
perturbation theory (golden rule). This corresponds to the classical picture of incoherent, 
sequential tunneling processes. On the other hand, there is experimental and theoretical 
evidence that in several regimes higher-order tunneling processes have to be taken into 
account. 

First, in the Coulomb blockade regime, sequential tunneling is exponentially suppressed. 
The leading contribution to the current is a second-order process in «o where electrons tunnel 
via a virtual state of the island. In Ref. [| the transition rate of this "inelastic cotunneling" 
process was evaluated at zero temperature. A divergence arises at finite temperature which 
requires a regularization. In Ref. |] this was done within an approximation which is valid far 
away from the resonances. In this regime, their results have been confirmed by experiments!. 

Second, at resonance, even though sequential tunneling occurs, higher-order processes 
have a significant effect on the gate-voltage dependent linear and nonlinear conduct anceHl. 
The energy gap between two adjacent charge states as well as the tunneling conductance is 
renormalized. Similar effects have been discussed for the equilibrium properties of the single- 
electron box A diagrammatic real-time technique developed for metallic islands™ as 
well as for quantum dots0§ allows a systematic description of the nonequilibrium tunneling 
processes. The effects of quantum fluctuations become observable either for strong tunneling 
«o ~ 1 or at low temperatures cto ln(E'c/T) ~ 1, where Ec denotes the charging energy (see 
below). The theory has been evaluated in Ref. §],0in the limit where only two adjacent charge 
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states are included (even virtually). Therefore, it was necessary to introduce a bandwidth 
cutoff ~ Eq. The predicted broadening of the conductance peak as well as the reduction 
of its height has been confirmed qualitatively in recent experiments on a single-electron 
transistor in the strong tunneling regime by Joyez et al@. However, a quantitative fit 
between theory and experiment requires introducing a renormalized value for the charging 
energy. 

In Ref. H we have used the real-time diagrammatic technique to obtain the current in 
second order in «o at low temperature when only two adjacent charge states are classically 
occupied, but all relevant virtual states are included. In this case no cutoff is required; all 
terms are regularized in a natural way. This analysis allows an unambiguous comparison 
with experiments where only bare system parameters enter. At resonance we obtain new 
contributions compared to the earlier theory of electron cotunneling. They emerge from a 
renormalization of the charge excitation energy and the tunneling conductance. For realistic 
parameters T/Eq ~ 0.05 and a}j = ~ 0.02 the corrections are of order 20%. 

In this article, we generalize the analysis of Ref. |^ in the following points: 

- We allow for arbitrary temperature. Our theory is, therefore, applicable both in the 
Coulomb blockade regime at low temperature as well as the classical regime at high tem- 
perature. While the effects of cotunneling are most dramatic in the first case, signatures of 
the higher order contributions are visible even in the second case. 

- We allow for arbitrary transport voltage, i.e., also extreme non-equilibrium situations are 
included in which more than two charge states are classically occupied. 

- We present analytic expressions for the current, and we discuss the asymptotic behavior 
at low and high temperatures. Renormalization group studies of the single-electron box0Jll 
predict a renormalization of the system parameters, which depends logarithmically on the 
relevant energy scale. This renormalization shows up in our formulas in new terms which 
are crucial at resonance. We are, therefore, able to extract the argument of the leading loga- 
rithmic terms. Finally, we compare with recent experiments^ and find very good agreement 
without any fitting parameter. 



This paper is organized as follows: 

In Sec. [Tl] we define the model, which is based on the standard tunneling Hamiltonian, 
and describe the diagrammatic technique which we use in the following. We, then, set up the 
general scheme for a systematic perturbation expansion in the dimensionless conductance 
a in Sec. fTJ. 

Sec. |IV| contains the general result for the expansion up to second order, in which we 
find four terms contributing to the total current 1^ = l[ 2 ^ + 1% + 1^ + I4 , Eqs. fl22f) - 



(Pq). The first one, J{ , is related to processes in which one electron enters and another one 
leaves the island coherently. This term is responsible for a finite transport in the Coulomb 
blockade regime. We recover the usual " cotunneling"0 but without any regularization prob- 
lems at finite temperature. The second and third contribution, I2 and I3 , describe the 
renormalization of the dimensionless conductance and of the energy gap for adjacent charge 
states, respectively. These terms are important at resonance. The fourth contribution, lf\ 
corresponds to processes in which two electrons enter or leave the island coherently. This 
plays a role in situations in which more than two charge states are classically occupied, i.e., 
at high temperature or high voltage. For a quantitative comparison with experiment all 
terms are important. 

The results simplify to closed analytic expressions in the two most interesting cases. 
The first is the low temperature regime (Sec. |V|) in which the renormalization of the system 
parameters plays an important role. The dimensionless conductance as well as the energy gap 
between adjacent charge states are effectively reduced. A discussion of this renormalization 
and a comparison is given in Sec. |VT| . 

The other case in which we obtain closed analytic expressions is the linear response 
regime for arbitrary temperature (Sec. |V11| ). We finally perform an expansion for low and 
high temperatures, respectively, and obtain very simple results for these limits in Sec. | V 1 1 1 . 

In general a numerical solution of a set of equations is necessary to evaluate the contri- 
bution I2 ■ 
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II. THE MODEL AND DIAGRAMMATIC TECHNIQUE 



The single-electron transistor (see Fig. |T]) is modeled by the Hamiltonian 

H = H L + H R + H 1 + H ch + H T = H + H T . (2) 

Here H r = J2kn e ln a lkn a rkn and Hi = J2qn e qncln c qn describe the noninteracting electrons in 
the two leads r=L,R and on the island. The index n = 1, ... N is the transverse channel index 
which includes the spin while the wave vectors k and q numerate the states of the electrons 
within one channel. In the following, we consider "wide" metallic junctions with N ^> 1. 
The Coulomb interaction of the electrons on the island is described by the capacitance 
model H ch = E c (h — n x ) 2 , where E c = e 2 /(2C) with C = C*l + Cr + C g . The excess 
particle number operator on the island is given by n. Furthermore, the 'external charge' 
en x = C^Vl + CrVr + C g V g , accounts for the applied gate and transport voltages. The 
charge transfer processes are described by the tunneling Hamiltonian 

= E E n>\knV l * + h.C. . (3) 
r=L,R kqn 

The matrix elements T^" - are considered independent of the states k and q and the channel n. 
They are related to the tunneling resistances i?x,r of the left and right junction by 1/Rt,t = 
(27re 2 /^)^ r (0)A^/(0)|T m | 2 , where iVi/ r (0) are the density of states of the island/leads at 
the Fermi level. The operator e ±l ^ keeps track of the change of the charge on the island by 
±e. 

We proceed using the diagrammatic technique developed in Refs. The nonequilib- 
rium time evolution of the charge degrees of freedom on the island is described by a density 
matrix, which we expand in Ht- The reservoirs are assumed to remain in thermal equi- 
librium (with electrochemical potential /z r ) and are traced out using Wick's theorem, such 
that the Fermion operators are contracted in pairs. For a large number of channels N, the 
"simple loop" configurations dominate where the two operators in a\ kn c qn from one term 
are contracted with the corresponding two operators c l La r fcn from another term H^, while 
the contribution of more complicated configurations are small. 
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The time evolution of the reduced density matrix in a basis of charge states is visualized 
in Fig. (Q) . The forward and backward propagator (on the Keldysh contour) are coupled by 
"tunneling lines" , representing tunneling in junction r. In Fourier space they are given by 

rl ; °exp[±/?(^-/i r )]-l 1 ; 

if the line is directed backward/forward with respect to the closed time path. Furthermore, 
we associate with each tunneling vertex at time t a factor exp (±iA n t) depending on the 
energy difference of the adjacent charge states 

A n = E ch {n + 1) - E ch {n) =E c [l + 2(n - n x )} . (5) 

If the vertex lies on the backward path it acquires a factor —1. We define a = J2r a o an d 

a(u) = £ r O! r (a;). 

The time evolution of the system is determined by an infinite sequence of transitions be- 
tween different charge states (see Fig. 0). This leads to a stationary probability distribution 
p n at the end. As shown in Refs. |5],|7| the probabilities p n are described by a master equation 
which reads in the stationary limit 

= [Pn'^n',n ~ Pn^n,n'] ■ (6) 

Here denotes the transition rate between n and n' . Only the latest transition, i.e., 

the rightmost irreducible part in the diagrams, is written explicitly. (We call a diagram 
part irreducible if there is no possible vertical cut separating the diagram into two parts.) 
The sequence of the preceding transitions is summarized in the stationary probability p n of 
charge state n. 

The tunneling current I T flowing into reservoir r can be expressed in terms of the corre- 
lation functions for the island charge 

C > (t,1?)=-i(e- i *®e i rt<)) i (7) 
C < (t,t')=i(e i ^e' i ^). (8) 
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Since the system is time-translational invariant, the correlation functions depend only on 
the time difference C(t,t') = C(t — £'), and we can define the Fourier transforms C(u) = 
J dte lujt C (t) . As shown in Refs. |6],|7| the tunneling current is 

/ r = -ie J du {o;+(a;)C > (a;) + a'^C^u)} . (9) 

For each contribution to the correlation functions the corresponding diagram contains two 
external vertices (representing the e ±llf terms in Eqs. (|7]) and (^) ) connected by a solid line. 
Two examples are shown in Fig. [5]. The value of the diagram is given by the probability of the 
states indicated on the left hand side of the diagram (p n in Fig. ||D multiplied with the value 
of the irreducible diagram part (consisting of two lines in Fig. ||]). The probability includes 
the sequence of preceding transitions which are disconnected from the latest irreducible part. 
While the calculation of the irreducible part follows straightforward from the diagram rules, 
we need some extra relation in order to determine the stationary probabilities p n . This 
relation can be the master equation Eq. (|6|). But instead of Eq. ([]) we use an alternative 
condition which is equivalent to the master equation. 

For this task we introduce the quantities C(u,n), I T (n) and I(n) = J2r^r{n) as the 
contribution to the correlation function and current in which the rightmost vertex describes 
a transition from n to n + 1 or vice versa. We can prove current conservation not only for 
the total current J2r Ir = but also for each part I(n) = 00, i.e., 

J du {a + (cj)C > (c;, n) + oT^C^uj, n)) = . (10) 

The normalization condition J2 n Pn = 1 together with Eq. fllCf ) determine the probabilities 
p n which can, then, be inserted in Eq. @ to calculate the current. 

III. SYSTEMATIC PERTURBATION EXPANSION 

For a systematic perturbation theory in powers of «o we write 

oo 

Pn= T,Pn ] 

k=0 
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(11) 



for the probabilities, and 

oo oo k 

C(u,n) =J2C {k) (u,n) = J2J2 c(l,k ~ l) (^,ri) (12) 

k=0 k=0 1=0 

for the correlation functions. The index k in pffi and denotes the term of the order 
O!o i n the expansion. The term C^ l ' k ~ 1 ^ is the contribution to the correlation function 
of the order k that is built up by a probability of the order I and an irreducible diagram 
part of the order k — I with < / < k. (One may write C"' k ~^ ~ p^a^iui) . . . a^iu^i) 
symbolically.) 

The conservation rule Eq. ( |T0D must hold in each order. This allows us to determine the 
probabilities recursively. If we know the probabilities of the order I < k — 1, then we get 
p^ from the k + 1-th order term of the normalization condition J2n Pn = $k,o an d Eq. ([II)]) . 

pW a +(A n )-pl k l ia ~(A n ) = 

i k 1 

— Y^fdu \a + (u;)C >{l > k - l \u, n) + a-(^)C <(/ ' fc - /) (cu ) ™)1 . (13) 
2m l=Q J l j 

Here, we have used 



C< k >°\uj,n) = 2mp<Zl 1 5{u- A n ) (14) 
C >(fe ' 0) (u;,n) = -27d#5(a;- A n ). (15) 

Having determined the probabilities up to pffi we get for the k + 1-th term of the current 
Eq. © 

/(*+!) = -ze | ^ {a+HC >(fc) (cu) + «;(cu)C <(fc) (cu)} . (16) 

In lowest order (sequential tunneling), we have k = in Eq. (|13|), i.e., p^ ^o; + (A n ) — 
pi°+iCt~(A n ) = with J2 n Pn^ = 1- The solution in equilibrium is = exp[— /^^(n)]/^ 
with Z = X)n exp[— pEch(n)]. At low temperature (and voltage) at most two charge states 
(n = 0, 1) are important, and all other states are exponentially suppressed. Higher-order 
processes modify the occupation, and the probabilities for the other charge states are only 
algebraically suppressed, but not exponentially. This leads to correction terms (n)^ k ' in the 
average occupation (n) = Y,n n Pn = J2k=o( n ) • 
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The current / = Jl = —In follows from Eq. ([16]) and is given in lowest order by 

' m = S K° ' + I/a(A») " /l(A„)] , (17) 

in which f T {uj) = l/{exp[/3(a> — /x r )] + 1} is the Fermi function of reservoir r. Higher-order 
terms account for new processes in which more electrons are involved as well as renormal- 
ization effects of the dimensionless conductance and of the energy gap of adjacent charge 
states. 

The rest of the paper is devoted to the calculation of the next order contribution (co- 
tunneling), i.e., k — 1 in Eqs. (|T3|) and flip]). We first calculate the correlation functions 
(u,n), insert the result into Eq. fll3"|) to get the corrections to the probabilities p^>. 
Then we have C^'°\u!,n) and CW(ai) from Eqs. (|HD, (|I5) and (|12D, and the current 1^ 
folloes from Eq. (Tl6|). 

IV. GENERAL RESULTS 

The first task is to evaluate all diagrams contributing to the correlation functions 
(u,n). For C < and C > there are 32 diagrams each. They contain one tunnel line and 
two external vertices. In Fig. |3] we show two such diagrams which contribute to C > ^°' 1 \uj, n). 
The other diagrams are generated from these two by putting the right and/or the left vertex 
of the tunnel line on the opposite propagator and/or change the direction of the tunneling 
line. In each case we have to adjust the charge states at the propagators in such a way that 
the rightmost (external) vertex describes a transition from n + 1 to n (or vice versa). By 
this procedure, we get in total 16 terms. The other 16 terms are obtained by mirroring the 
old diagrams with respect to a horizontal line and changing the direction of all lines. The 
diagrammatic rules imply that the mirrored diagrams are minus the complex conjugate of 
the old ones. 

The first and second diagram in Fig. |3| have the values 

pf / dJ a + (u') \ - 1 (18) 
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and 



/ dw' a + (u/) = . . (19) 

n J V J u-u' + iO+ (uj- A n + i0+) 2 V ; 

The integrands have poles. But these poles are regularized in a natural way (this follows from 
our theory and has not to be added by hand) as Cauchy's principal values and delta functions, 
l/{x + iO + ) = P[l/x]-iir6(x), and their derivative l/{x + i0 + ) 2 = -(d/dx)P[l/x]+iir6'(x). 
The integrals are well-defined at these points. Another condition for the existence of the 
integrals concerns the behaviour at \uj'\ — > oo. Since a(cu') raises linearly at large energies, the 
integrands of some integrals may not decay fast enough at \uj'\ — > oo to ensure convergence 
for a single diagram. The sum of all diagrams, however, is convergent. To see this, we 
use for each diagram a Lorentzian cutoff U 2 /{U 2 + u' 2 ) with U larger than all relevant 
energy scales. Then, all intermediate results become convergent but may be ^/-dependent. 
To be more precisely, the cutoff enters only in the form W = ln(|^. While the W is 
present at some intermediate steps (even in the probabilities p^), we find that the physical 
quantities like the average charge and the current are independent of W. This means that 
the divergencies are only an artefact of considering single diagrams and are not present for 
the sum of all contributions. 

Defining 0„ ir = «o(A„ - // r )Re\I> («Jp(A n - /i r )) and n = E r 0n,r, where * denotes 
the digamma function, and making use of A n+1 — A n = 2E C , we obtain the lengthy but 
complete result 

C >( -°' r >(u,n) = 2m 

p ( °W(A n + A n _i - u) +pfa-(u) 



Rc 



+ 



p^a + (A n+1 + A n -uj) + pf +1 a-{u)] Re 
+ P i°2 1 5(c-A n _ 1 )^ 

-%-A n )[piV0n-l+P^ 



-R n (u) 



u - A n + i0 + 



(0) x , A \ 0n+l — ' 
— T Pn+1^ ~ V) ^ 



+ 



'n+1 



)+P { nlid<f> n +p^5\uj-A n )\2<p n 



>n-l\ 



+a wA [S(u - A n ) - 5(cu - A n _x)] + a Wp% [6(u - A n ) - 5(cu - A n+1 )]} . (20) 
Here d stands for 
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'Trh~~J^~ and we have used the abbreviation 

2Ec on x 



1 



1 



R n {uj) = — 

v ; uj-A n + iO+ w-A n _i+iO+ 



(21) 



The other contribution to C > ^(cu, n) = C > ^°' 1 \uj, n) + C > ^ 1 '°\uj, n) involves probabilities 
in higher order, see Eq. ([15]) for k = 1. 

The expression for C < ^(cu, n) looks very similar. It is obtained from C > ^(to, n) by 
changing the direction of all lines and adjusting the charge states. This implies a + <-> a~, 

A„-l <-> A n+ i, <-> n+ i, 2£ C «-> -2-^C, Pn-l ^ Pn+2, Pn <-> Pn+1, #n <"> ~Rn+l, and, 

furthermore, there is a global minus sign. 

We can solve Eq. ([13|) with the help of C >< -°' 1 ' ) (c(j, n) and C^ ' 1 '^, n) for the proba- 
bilities pffl. Due to the last line of Eq. (p0|) the corrections of the probabilities, pW, are 



,(°) 



,(°) 



As a 



VF-dependent. But the VF-dependence is of the form — a^W 2pffl — p n +i — Vn~\ 
consequence, in C^>{oS) = Y^ n ^ { w i n ) as we ^ as m (^r the ly-independent terms drop 
out, i.e., there is no divergency problem for the physical quantities. 

Now we are ready to formulate the most important result of this paper. Inserting the 
result for C^(uS) in Eq. flTop with k — 1 we obtain the cotunneling current which can 
be divided into four parts, 1^ = Y^t=i h with 



R n {uf + R n+1 {uof 



A (2) =Ep1 0) / dujB(u)a(u)ReR n (u) 2 , 

L n J 

n 

4 2) = -\ E {pf +pSi) / dua(u)Re [R n (u) - R n+1 (u)} , 



(22) 

(23) 
(24) 



r(2) 



47r 2 e 



^L + M"fi( A n) D 

u;J — r—r Re 



q L (u)a R (A r 
a(A n ) 



-Re 



- A„ + i0+ 
1 



RJu) 



+p^a + (A n+1 + A n - w) tt—, Re 



~Pn+2 a (A„+l + A n - CJ 

- (L <-> R) . 



a(A n ) 
a(A n ) 



Re 



— A„ + z0~ 
1 

w - A n + i0+ 
1 



u; - A n + iO 



^Rn(u) 
R n+ i(uj 
-R n+ i(uj) 



(25) 
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Here we used the definitions 

_ A7r 2 ea L (uj)a R (uj) 
B{u) = — -j- [f R {u) - f L (u)\ (26) 

and 

p n = p<? - a W [2pW - pffij - A] - (p?) - pLt)^ + (pffij - pf )<90„ . (27) 
In the low temperature limit, B(Ao) is the sequential tunneling result. It enters in l[ 2 \ 

(2) 

and J3 in different ways. 

The first term /) describes processes in which one electron enters and other one leaves 
the island coherently. In the Coulomb blockade regime only one charge state is occupied 
in lowest order at low temperature, p^ = 5 n>0 . At T = 0, the integrand is zero at the 
poles of R n (u), and we can omit the term i0 + . In this case, we recover from Eq. fl2~2~D 
the usual "cotunneling" result^. At finite temperature, however, the term +i0 + is needed 
for regularization which is not provided by previous theories^. The result Eq. ( p2"|) is also 
well-defined for T^O. 

(2) (2) 

Furthermore, we are able to describe the system at resonance. In this regime, Ir, and 1% 
become important. These terms describe the renormalization of the tunneling conductance 
«o and the energy gap A , respectively. They can be traced back to vertex correction and 
propagator renormalization diagrams (see Fig. 0). We will discuss these contributions in 
more detail in Sec.|V| and Sec.[V]. 



At temperature or bias voltage, there are processes in which two electrons enter or leave 
the island coherently. After such a process, the island charge has increased or decreased by 
2e. These processes are described by the fourth term I\ , Eq. (p5|). 

All integrals in Eqs. ( p2"D - (|25| ) can be performed analytically. We do not write down the 
lengthy results for the general case here. They can be found for the most interesting limits 
in Sec. [V] and |V11| . The quantities p n which enter Eq. (|23| ) are determined by J2 n Vn = and 

p n a + (A n ) -p n+ ioT(A n ) = 
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r,(°) _ p /3(A„+A„-i-/ir-^)J°) 

E Pn-l e Pn+1 r 
e /3(A n +A„_i-/V-/v) -1 a ° 

r,r' 

A„-i - A*r , , A„, - /i 



x (-9 [(A n - /i r )0„_i, r ' + (A n _i - /ir.)0„. r '] ^—^ — —4> n -iy + , . 



T 4>n,r' 



7,(0) _ p^+i+A^-^-^,) (0) 

+ 2^ — — — — : " 



e /3(A„ + i+A n -/.i r -/i r /) _ i 

A n — fi r A n+ i — fl r 
; <?>n,r' H ^T - 

#<9a + (A„) - p^ 1 9a-(A n )l [20 n - (j) n+1 - </> n _i] (28) 



X {—B [(A n+1 — ll r )(j) n y + (A n — ft r )4> n +iy} " ^ <fin,r> H — </>n+l,r' 



In principle, the problem is now solved. One can solve Eq. ( p8|) for p n numerically and 
will then get the average charge and the current. In Figs. (| and [| we show the differential 
conductance as a function of the transport voltage for two different values of the gate charge. 
There are steps at voltages at which new charge states become important. Due to cotun- 
neling the steps are washed out. Furthermore, there is a finite conductance in the Coulomb 
blockade regime although the sequential tunneling contribution is exponentially suppressed. 

But we can do even better and find analytic expressions for the two most interesting 
cases. The first one covers the regime in which only two charge states are occupied in lowest 
order (then four charge states have corrections in first order). In this regime the Coulomb 
oscillations are most pronounced and the system is expected to behave like a multichannel 
Kondo model. But to study the crossover from the Coulomb blockade to the classical (high 
temperature) regime, one has to allow for all temperatures. This is done in the second case, 
in which we consider the linear response. 



V. LOW TEMPERATURE NONEQUILIBRIUM TRANSPORT 

At low temperature, ksT <C Eq, and a transport voltage which may be finite but small 
enough, only two charge states (n = 0, 1) are occupied in lowest order, p^ ps for n ^ 0, 1. 

Sequential tunneling gives for the average charge just (n)^ = pf = a + (A )/a(A ), and 
the current is I^ 1 ' = B(A ). 
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The probabilities in next order, p£\ are determined by Eq. (0). In the low temperature 
limit, the first and second term in Eq. ( p8|) are exponentially small. We find 



p9 



(o). 



0-1 - 0l] + (Pi 



(0) 



(pi 0) - j9 o o) )«90 o + p^dfa + dpf [20 o 



(0); 



,(o) 



- 2$ ) )a W 
2pf ) )a W 



P2 



(1) 



p^dfa + p^'a W 



(o). 



(29) 
(30) 
(31) 
(32) 



while all other corrections vanish, and therefore 



(n) 



(i) 



1 d 



(<P-i-M+P [ i" (00 -<k) 



(0) 



2E C dn K 

The total second order, "cotunneling" current Eqs. (|22"D - 



(33) 



simplifies to 



n 



(2) 



duo B(uo)a(uj) Re 



pf Roico) 2 + pT> R 1 (u 



(o) 



r(2) 



-^B(Aq) I duoa(uj)Re R (uj) 2 + R^ 2 



r(2) 



and /I 2 '* 

r(2) 



1 &B(A 
"2 <9A 

0. All integrals in Eqs. 

47r 2 e 



(34) 
(35) 
(36) 



du a(uj) Re [R (uj) — Ri(u)] , 

can be performed analytically. The result is 



h 



P 



(0) 



A„_l-/iR A n -/iR 

-0n-l,L (pn.L 



Er 



n=0,l ^ E C 

+8 [(A„_x - // R )0„_l, L + (A n - // R )0 n ,L] 



(L «-> R) 



(2) 



5( A ) 1 ^ 01 + <9[20 o + 0_ x + X ] 



I 



(2) 



dB(A c 
dA 



(20o - 0_i - 00 



(37) 
(38) 

(39) 



I^ XJ = 0. The only contribution is the term with n = in Eq. ([34"1). This gives 



In the Coulomb blockade regime, we have p^ = 1, p-f 1 = and 7^(Aq) = 0, which 
yields 4 1} - r(1) 

the well-known result of inelastic cotunnelin 

At resonance (p^ ^ p^), the terms describing the renormalization of the tunneling 

(2) (2) 

conductance «o and the energy gap A , an d -^3 ■> become important. Suppose that the 
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current can be expressed by the sequential tunneling result with renormalized parameters a 
and A plus some regular terms, I(a , A ) = I^(a, A) + regular terms. The expansion up 
to 0{al) of 

/«(*, A) = ^/ (1) («o, A ) + (A - A )%M . (40) 
a dA 

and the identification with Eqs. (|38|) and fl39f) yield the lowest-order corrections for the 
renormalized quantities. 

In Figs. ^ and |7| we show the second-order contribution to the differential conductance 
G = dl/dV for the linear and nonlinear regime (in the following we choose = ct^). The 
dimensionless conductance a as well as energy gap A is decreased in comparison to the 

(2) 

unrenormalized values «o and A . For this reason, is negative since the conductance is 

(2) 

reduced, and 1% is positive since the system is effectively "closer" to the resonance. 

In Figs. |8] and || a comparison of the first order, the sum of the first and second order, and 
the resonant tunneling approximation^ (where the cutoff is adjusted at Eq) is displayed. 
The deviation from sequential tunneling is significant and of the order 20%. The agreement 
with the resonant tunneling approximation provides a clear criterium for the choice of the 
bandwidth cutoff. We have checked the significance of third order terms ~ ajj by using the 
resonant tunneling formula&B and exact results for the average charge in third-order at zero 
temperature^!. For the parameter sets used in the figures, the deviations to the sum of first 
and second order terms were smaller than about 2%. Therefore, at not too low temperatures, 
second-order perturbation theory is a good approximation even if the tunneling resistance 
approaches the quantum resistance. 



VI. COMPARISON WITH RENORMALIZATION GROUP RESULTS 

Renormalization group techniques have been applied for the single-electron box (which 
is equivalent to the transistor at zero bias voltage) including two charge statesS0. Starting 
with some initial high-energy cutoff (which is of the order of the charging energy Eq), 
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one integrates out the high-energy modes. This leads to a renormalization of the system 
parameters. This procedure has to be performed until the cutoff reaches the largest energy 
scale toe which is important for the lowest order tunneling processes. This leads to 

A = A = 1 (41) 

a A l + 2a \n(E c /uJc) 

Using the resonant tunneling approximation!'! we generalized this to non-equilibrium situ- 
ations and find (for Eq larger than any other energy scale) 

a 



a 



A 



l + 2E r «& [ln(fk)-Retf (? 



2tt 



A + 2 E, <*o li, 


Mt) 


Re*(z^=^)] 


l + 2Er<*5 




- 







(42) 
(43) 



But in both cases, the result is rather qualitatively, and the constants in the arguments of 
the logarithms remain unknown. Here, in the cotunneling theory, we have taken all relevant 
charge states into account, i.e., the cutoff is provided naturally by the energy of the adjacent 
charge states and has not to be introduced by hand. Therefore, we can extract the complete 
arguments from Eqs. fl38|) and (j39|), 



a 



On 



■1 + ln 



7T 



9 ((A \p , T , ^(Ao-^) ' 



A = A -2]Ta^(Ao-/i r 



2tt 



(44) 
(45) 



In the following we discuss different limits. 



A. Linear Response 



At zero bias voltage the system is "on resonance" if (3Ao <C 1 and "off resonance" if 
/3A > 1. 

While the resonant tunneling formula yields "on resonance" 



a A 



«o A 1 + 2a [7 + In (^) 



(46) 
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we get from cotunneling 

a = ao 
A = A 



1 - 2a y-l + 7 + In 
1 - 2a ( 1 + 7 + In ( 



'(3E< 



c 



V 71 

13 E c 



IT 



(47) 
(48) 



We see that the temperature is the relevant energy scale which determines the renormaliza- 
tion. 

In the opposite case, i.e., "off resonance" the resonant tunneling approximation leads to 

A 1 



a 



a 



A l + 2a ln(f)' 



(49) 



and the cotunneling yields 



o. = Qq 
A = A 



1 - 2a i^-2 + In 
1 - 2a (l + In 



2E, 



A 
2Ec 

A 



(50) 
(51) 



Now, the level splitting itself provides the energy at which the renormalization procedure 
stops. 

The conductance as well as the energy gap is reduced during the renormalization, i.e., 
the system is effectively "closer" at the resonance but with higher tunnel barriers than the 
bare system. As a consequence, is negative and is positive. 



B. Nonlinear Response 

At finite bias voltage a new energy scale comes into play. Furthermore, since the Fermi 
levels of the reservoirs are separated, the relevant energy scale may depend on which reservoir 
is responsible for the renormalization. For /% — /xl = eV ^> ksT, being on resonance with 
the left Fermi level means that k B T ^> |A — /j L |, i.e., that again the temperature provides 
the dominant energy scale for the left reservoir, but |A — /Jr| ~ eV is the relevant energy 
for the right one. 

The resonant tunneling formula leads to 
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a 



a 



l + 2a L [ 7 + ln(^)]+2^ In ) 



A-/x L 

and the cotunneling gives 



A - fj, L + 2a$eV In 


(30 




1 + 2^ 


[7 + ln(^) 


; + 2a* In (f ) 



(52) 
(53) 



lot 



7 + In 



'PE C 



71 



2c# I -2 



In 



/2£r 



V eV 



A - /x L = A - /x L - 2«o (A - /x L ) 



1 + 7 + In 



7T 



1+ln 



/2R 



V eV 



(54) 
(55) 



In the opposite case, i.e., |A — /x L | k B T and |A — /XrJ 3> /cb^ the resonant tunneling 
yields 

a 



A + 2E r ^/x r ln(^) 



(56) 
(57) 



in comparison to cotunneling 

a = «o 



A = A -2^^(A -/x r 



1 + ln 



2R 



(58) 
(59) 



|A - /x r |, 

The renormalization of the energy gap differs qualitatively from that in the linear re- 
sponse regime, since now there is a competition between the renormalization of the system 



towards the left and the right Fermi level. This effect is displayed in Fig. [T(] which shows 
the differential conductance as a function of the external gate charge. The position at which 
the system is in resonance with the left Fermi level is shifted towards the right Fermi level 
with increasing coupling strength to the right reservoir. 



VII. LINEAR RESPONSE FOR ARBITRARY TEMPERATURE 



In the last section we discussed the effect of the quantum fluctuations at low temperature, 
when Coulomb oscillations are pronounced. In this case, at most two charge states are 
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occupied classically. But also at high temperature, the asymptotic behavior shows deviations 
from the sequential tunneling result. Furthermore, for a comparison with experiment, it is 
necessary to describe also the crossover between the low and the high temperature regime. 
In this section we discuss the linear response, i.e., we set V — 0, but we allow for all charge 
states. 

The corrections to the probabilities are then 



P 



P = (pf - A) dK-X - (Pffil - 9<Pn ~ (0.-1 - 0r. 



(0) 



V-l — H>n 



)-a w(2p^- P ^ 1 - P ^ 1 



(60) 



which implies 



(n) 



(i> 



J d_ 



Y.P 



(o) 



>n-\ — Vn) 



(61) 



A systematic perturbative expansion of the partition function (up to order aft) was performed 
by Grabert^l. The result Eq. ( |6"ID is identical to his result in order «o, which at T = reads 
(n>« =a ln[(l + 2n x )/(l-2n x )]. 

Using the definition S(A n ) = dB(A n )/dV we find the differential conductance = 
dl^/dV = Y?i=iGf ] with 



J n 

G? ] = ~E S ( A n) / dwa(«)He{(pW + [^(a;) 2 + it! n+1 (a;) 2 " 

^ n ^ 

H-rf^M + f^ x R^M - (3 (pf + pi ],) Ep^'M] 

n' J 

= -i E (p1 0) + p^) / ^ «M Re f^H - 



(62) 



(63) 
(64) 
(65) 



In comparison to the previous section, a new term, G\ , appears, which is related to new 
processes where the island charge is increased or decreased by 2e. Again, all integrals Eq. (| 
- (6"5p can be performed analytically. While the sequential tunneling result gives 
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G« 



(0) , (0) 



/?A n /2 
sinh j3A n 



(66) 



we find for the cotunneling 



(2) 



Gas 



A 



n-\<Jipn-l 



A 



Er 



n^Yn Yn 



>n-X 



Er 



(67) 



G a 



E 



sinh /3A r 



D (o) , „(o) 



->n+l 



+ <9(20„ + 0„-i + n +l) 



G 



(2) 



E 



d / f3A n /2 
dA n Uinh/?A r 



G 



(2) 



G a 



^E- 



Pn-1 + Pn+X 



, A n _x + A r 



sinh/3 (A„_i + A n ) 



-A n _i90 n - A„( 



2£ f 



(68) 
(69) 
(70) 



In Fig. |TT| we compare our results with recent experiment s@. The temperature depen- 
dence of the Coulomb oscillations were measured for a sample with a conductance ao = 0.063. 
Our results in second-order perturbation theory agree perfectly in the whole temperature 
range. For stronger tunneling higher-order effects such as the inelastic resonant tunnelingB 
would be relevant. 

We emphasize, that only bare values for a and Eq have been used here, as determined 
directly in the experiment. In contrast, the resonant tunneling approximation with bare 
values of the charging energy would lead to a deviation from the experiment by about 10%. 
Thus, the inclusion of higher-order charge states within second-order perturbation theory, 
as presented in this paper, is important for comparison with experiments. 

While at resonance the new terms are crucial, the Coulomb blockade regime is sufficiently 
described by Eq. ( |3~4D which yields a good agreement between theory and experiment. 



VIII. ASYMPTOTIC BEHAVIOUR 

In linear response, the whole temperature regime is covered by Eqs. ( |66"D - ([701). At high 
and low temperatures these formulas become very simple. 
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A. High-temperature expansion 



At high temperature (flEc <C 1), Coulomb oscillations are washed out, i.e., there is no 
gate charge dependence of the conductance, and using Euler-MacLaurin formula we can 
replace the sum in Eq. by an integral. In this way we average the conductance over 
all possible gate charges, loosing the information about the difference between minimal and 
maximal conductance. In linear response we have = exp[— PEc(n — n x ) 2 ]/Z. We expand 
(/3A n )/ sinh(/?A n ) in (3 A and perform the integral which approximates Eq. This yields 
G^/Gas = 1 - (3E c /3 + (PE C ) 2 /15 + 0((PE c ) 3 ). 

For the cotunneling contribution we expand the digamma functions which enter n = 
a A n Re* (i^) using Re^(ix) = Re^(l + ix) = -7 + C(3)x 2 - C(5)x 4 + C(x 6 ), and we 

((3E c ) 2 ((3)/n 2 and Gf = 0, i.e., 

G . fiEc 



find G\ 



(2) 



G? 



G a 



1 - 



+ ((3E C ) 2 



1 ^ 6C(3) 
77 + — — a o 

15 7l Z 



+ 



En? 



(71) 



This result was also derived in Refs. 19.20 



B. Low-temperature expansion 

At low temperature {(3Eq ^> 1), Coulomb oscillations appear. The maximal conduc- 
tance is reached at half integer values of n x . In sequential tunneling the maximal con- 
ductance saturates at one half of the asymptotic high temperature value. The quantum 
fluctuations, however, lead to a further reduction. Approximating the digamma function by 



Re W(ix) = In \x\ for x ^ 0, we find G^J G as = -a , G^Lx/ G ^ = «o 1 - 7 - In 



and Gr max = GT> = 0, i.e., in total 



Gv 



G aK 2 



77 - «o 



7 + In 



'0E, 



c 



7T 



+ O (a 2 ) (72) 



with 7 being Euler's constant and G as = 1/(Rl + Rr) the asymptotic high temperature 
limit. The peak conductance depends logarithmically on temperature. This result may 
be interpreted as a renormalization of G as or ao&~&0El. It shows a typical logarithmic 
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temperature dependence since, at least in the equilibrium situation, the low-energy behavior 
of the system is expected to be that of the multichannel Kondo mo deM 

In Fig. [12] we compare the limiting expressions Eq. ([H]) and (|72|) with the full second- 
order result Eq. fl66|) - flT0"|). Our theory interpolates perfectly between the logarithmic 
behavior at low temperature and the high-temperature result. 



IX. CONCLUSION 

In conclusion we have evaluated the current of the single-electron transistor consistently 
up to second-order perturbation theory (cotunneling) at arbitrary temperature and trans- 
port voltage. The approach is free of divergences and provides cutoff-independent results. 
Out of resonance we recover the usual "cotunneling" contribution at low temperature!. At 
resonance, however, we find new terms which are significant for experimentally realistic pa- 
rameters. They describe the renormalization of the dimensionless conductance and of the 
energy gap of two adjacent charge states. Both parameters are reduced in comparison to 
the unrenormalized value and depend logarithmically on temperature and bias voltage. We 
can extract the arguments of the logarithms which are not provided by renormalization 
group techniques and which are important for comparison with experiments. Furthermore, 
we include processes in which two electrons enter or leave the island coherently. The corre- 
sponding term is important to describe the deviation of the sequential tunneling result at 
high temperatures. 

We have derived analytical expressions at low temperatures, including nonequilibrium 
effects. In addition we have studied the linear response regime for arbitrary temperature. We 
are, thus, able to describe the crossover between the Coulomb blockade and classical high- 
temperature regime. A comparison with experiments shows good quantitative agreement. 

We like to thank D. Esteve, H. Grabert, and P. Joyez for stimulating and useful discus- 
sions. Our work was supported by the "Deutsche Forschungsgemeinschaft" as part of "SFB 
195". 
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FIGURES 
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FIG. 1. Single-electron transistor. 
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FIG. 2. A diagram showing contributions to sequential tunneling (Sqi an< ^ ^io) an< ^ c °t un - 

(2) 

neling: usual cotunneling with two electrons tunneling coherently (contribution to £qo)' vertex 

(2) (2) 

correction (Sq {), and propagator renormalization (£•[ q). 
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FIG. 3. Two diagrams contributing to C > ^°' 1 \u;, n). 
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FIG. 4. The differential conductance as a function of transport voltage for T/Ec = 0.05, 
n x = 0, and «o = 0.01: sequential tunneling and sequential plus cotunneling contribution. 
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FIG. 5. The differential conductance as a function of transport voltage for T/Ec = 0.1, 
n x = 0.5, and uq = 0.01: sequential tunneling and sequential plus cotunneling contribution. 
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FIG. 7. The second-order contribution to the differential conductance = Y$=i for 
T/E c = 0.01, a = 0.02 and V/E c = 0.2. 
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FIG. 8. The differential conductance for T/Eq = 0.05, ao = 0.04 and V = 0: sequential 
tunneling, sequential plus cotunneling, and resonant tunneling approximation. 
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FIG. 9. The differential conductance for T/E c = 0.01, a = 0.02 and V/E c = 0.2: sequential 
tunneling, sequential plus cotunneling, and resonant tunneling approximation. 



28 



0.0 

0.40 0.45 0.50 0.55 0.60 

FIG. 10. The differential conductance T/E c = 0.01, V/E c = 0.2, aft = 0.0001, a% = 0.01 
(solid line) and = 0.001 (dashed line). 




T/E c 

FIG. 11. Maximal and minimal linear conductance for Eq = IK and ao = 0.063. The dashed 
line is sequential tunneling, the solid line sequential plus cotunneling. The data points are experi- 
mental data from Ref. 
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FIG. 12. Maximal linear conductance for Eq = IK and ao = 0.063. The solid line (full first 
plus second order result) interpolates between the long-dashed line (low-temperature expansion) 
and the dashed line (high-temperature expansion). The data points are experimental data from 
Ref. |i~6. 
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